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Abstract 

The coincidence problem for planar patterns with A^-fold symmetry is 
considered. For the A^-fold symmetric module with < 46, all isometrics of 
the plane are classified that result in coincidences of finite index. This is done 
by reformulating the problem in terms of algebraic number fields and using 
prime factorization. The more complicated case > 46 is briefly discussed 
and = 46 is described explicitly. 

The results of the coincidence problem also solve the problem of colour 
lattices in two dimensions and its natural generalization to colour modules. 



1 Introduction 

The concept of coincidence site lattices (CSLs) arises in the crystallography of grain 
and twin boundaries [T . Different domains of a crystal do have a relationship: There 
is a sublattice common to both domains across a boundary, and this is the CSL. 
This can be seen as the intersection of a perfect lattice with a rotated copy of it 
where the set of points common to both forms a sublattice of finite index, the CSL. 
Up to now, CSLs have been investigated only for special cases, for example for cubic 
or hexagonal crystals |2] . With the advent of quasicrystals infinitely many new cases 
arise: quasicrystals also have grain boundaries, and one should know the coincidence 
site quasilattices OH]. In a rather different context, multiple coincidences of families 
of ID quasicrystals have been applied in constructing quasicrystals with arbitrary 
symmetry (in higher dimensions) jSlEI- An application of these results was made 
by Rivier and Lawrence to crystalline grain boundaries, which themselves turn 
out to be quasicrystalline. This is an important example of the relevance of a 
coincidence quasilattice. The experimental evidence was provided indirectly by Sass, 
Tan and Balluffi in the 1970s jSj, but beautifully by the observations of growth of 



quasicrystalline grain at the grain boundary between two crystals by Cassada, Shiflet 
and Poon |9| and by Sidhom and Portier JU]. Gratias and Thalal jTT], on the other 
hand, used quasicrystal concepts in a different context to embed the two crystal 
grains adjacent to a grain boundary in a higher dimensional perfect lattice. So an 
extension of the CSL analysis to more general discrete structures is desirable. 

In this paper we give a unified treatment of the coincidence problem for planar 
structures with general A^-fold rotation symmetry, extending previous [31 and 
parallel p2] work and putting it in a more general setting. This is what is needed for 
quasicrystalline T-phases which are quasiperiodic in a plane and periodically stacked 
in the third dimension. Icosahedral symmetry in 3D requires different methods and 
will be described separately j^. Common to both is the necessity of an attack in 
two stages: not only do we have to find the coincidence isometries (the universal part 
of the problem), but also the specific modifications of the atomic surfaces (also called 
windows or acceptance domains) that are needed to describe the set of coinciding 
points. 

In order to describe this scenario, we start with the coincidence problem of the 
square lattice 1? . The set of coincidence transformations for 1? forms a group, the 
generators of which can be given explicitly through their connection with Gaussian 
integers. Simultaneously, the so-called S-factor or coincidence index can be calcu- 
lated for an arbitrary CSL isometry. Though this is not new, the approach we use 
here can be generalized to quasiperiodic planar patterns with iV-fold symmetry. The 
description of this more general case and the tools necessary to tackle it is the main 
aim of this article. 

In two dimensions the classification of CSLs is the same as the classification of 
colour lattices with rotational symmetry ([TSj, Section 5.8). In that setting the 
S-factor or coincidence index is the number of colours and the different coloured 
sublattices are the different cosets of the CSL in the original lattice. This is because, 
as long as the symmetry group consists of 2D rotations, all members of the symmetry 
group commute with the CSL rotation, thus ensuring that the CSL is invariant under 
the symmetry group. For indecomposable groups in higher dimensions no non-trivial 
orthogonal transformation commutes with all symmetries so there is no longer this 
equivalence. The only non-trivial rotation groups of 2D lattices are C3, C4 and Cg. 
The prime numbers p for which there exist p-colour lattices with these symmetries 
are listed in [T3], p. 76, and coincide with the sets of primes in the denominators 
of the Dirichlet series given at the end of Section 0] for the cases n = 3 and n = 4. 
(Note that N = 2n for n odd and N = n otherwise as will be explained later.) For 
non-lattices the solution of the CSL problem in 2D can be regarded as a classification 
of colour modules in the plane. An r -colour n-module is a pair of n-modules (A^, 
M-i) such that M.I has index r m M. and is invariant under the symmetry group of 
M.^ while no other coset of A^i has this property (see Section 3 for a definition of 
n- modules). The colour of a point in M. is then determined by its coset mod M.i. 
In this light, the results of Sections El and El can be interpreted as finding, for each 
n, the numbers r for which there are r-colour n-modules and what these n-modules 
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are. 

The paper is organized as follows. In Section 2, we review the coincidence problem 
for the square lattice and formulate it in terms of Gaussian integers. This enables 
us to describe the group structure and the coincidence indices explicitly and to 
introduce the concepts needed for the generalization in Section 3. There, the main 
structure is derived with the aid of the algebraic number theory of cyclotomic fields, 
followed by various explicitly worked out cases in Section 4. They include 8-, 10- 
and 12-fold symmetry, the most important cases for quasicrystalline T-phases, and 
thus cover all cases linked to quadratic irrationalities [16> . In Section 5 we then 
show, in an illustrative way, how to use the method for the eightfold symmetric 
Ammann-Beenker rhombus pattern and the tenfold symmetric Tiibingen triangle 
tiling. We give an explicit formula for the necessary correction of the coincidence 
index. In Section 6 we discuss certain details to be dealt with for > 46, where 
the variety of modules rapidly increases, though this does not affect the generality 
of our findings. The case = 46 (n = 23) is presented in some detail. This 
is followed by some concluding remarks, while the two appendices cover further 
examples (Appendix A) and proofs of technical results used in Sections 3 and 4 
(Appendix B). 



2 The square lattice: a warm-up exercise 

Let us consider the CSL problem for the square lattice 1? . We focus on pure rota- 
tions first and deal with the easy extension to reflections later. Consider therefore a 
rotation (i.e., an element of the group 5*0(2) = 50(2, M)) and ask for the condition 
that it maps some lattice point to another one. Clearly, rotations through multiples 
of 7r/2 do this. They form the cyclic group O4 — an index 2 subgroup of -D4, the 
point group of 1? . 

But there are more cases, as can already be seen from the growing number of 
lattice points on expanding circles, summarized in the coefficients of the theta- 
function of the lattice, cf. (T7j, 



= l+^r(M)a;^ (1) 

M=l 

= l + Ax + Ax^ + Ax'^ + 8x^ + ... 

Here, '&3{x) = J2qez -^'^'^ = 1 + 2x + 2x'^ + 2x^ + ... is Jacobi's theta-function and 
r(M) denotes the number of integral solutions of the equation + 6^ = M, see jTHj 
for details on r(M). This number is only slowly increasing but is unbounded, so 
there is an infinite number of rotations that map one lattice point to another. 

As is obvious (cf. ^H] and references therein), the set of coincidence rotations 
(or CSL rotations) consists of all rotations R through angles (f with sin{(f) = a/m 
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and cos{(f) = b/m rational, and hence is identical with the group S0{2,Q). This 
requires integral solutions of the Diophantine equation 

a^ + b^ =m^ , (2) 

where we need consider only the primitive solutions, i.e., gcd(a, 6) = 1. They are, 
of course, given by the primitive Pythagorean triples [^. For a primitive solution, 
the set of coinciding points forms a sublattice of of index m, whence 1/m is the 
fraction of lattice points coinciding. We call m the coincidence index of R, denoted 
by S^2(-R), or for short. This index is often called the S-factor |2llSlll]. 



The number of CSL rotations with given index 

Without determining the rotations explicitly we can calculate their possible indices 
and the number of different rotations with each index as follows. 

The number of primitive solutions of Eq. (j2)) can be derived from the well-known 
formula (cf. ^) 

r{M) = A{di{M) - d^{M)) , (3) 

(where dk{M) counts the number of divisors of M of the form 4i + k) for the total 
number of integer solutions of 

a^ + b^ = M . (4) 

If we write M = 2^MiM^, where Mi and M3 are maximal divisors of M composed 
of primes congruent to 1 or 3 (mod 4), respectively, then Eq. (jH)) can be equivalently 
expressed as 

f4rf(Mi), if M3 is a square, 
r[M) = < (5) 
1^0, otherwise, 

where d{Mi) counts all the divisors of Mi. When (as in our case) M is a square, 
the first alternative in Eq. occurs. The number of primitive solutions, r*(m^), 
of Eq. can now be derived from the "input-output" principle (cf. [TH], Thm 260) 
as 

where p runs through all prime factors of m, pp' through all pairs of distinct prime 
factors of m, and so on. After substituting Eq. (jSJ in the right hand side of Eq. (jH)) 
and then counting the contributions of the factors of m one at a time, it can be seen 
that 

2n I 4(i*(m), if m has prime factors = 1 (4) only, 
r*{m ) = < (7) 
[0, otherwise, 

where d*{m) counts the squarefree divisors of m. We note that the number of CSL's 
(as distinct from CSL rotations) of index m in the square lattice is a quarter this 
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number, since each is itself a square lattice stabilized by the rotation group of the 
square (of order 4). (Note however that not every sublattice with square symmetry 
is a CSL.) 

So far we have: 

Theorem 1 The coincidence indices of the square lattice are precisely the numbers 
m with prime factors = 1 (mod 4) only. The number of coincidence rotations f{m) 
with a given index m is 

/(m)=4d*(m) (8) 
and the number of CSL's with index m is 

f{m) = d*{m). (9) 



CSL rotations and Gaussian integers 

We have settled the question of what numbers occur as coincidence indices of CSL 
rotations of and how many rotations there are with each index, but there is still 
more to be said. 

We have seen that the set of CSL rotations forms a group (5*0(2, Q), in fact). 
Let us introduce the notation 

SOC{Z^) := {R G 50(2) I S(i?) < oo} (10) 

for it. We shall investigate its structure and derive independent generators. 

The most transparent proof of Eq. Q (that given in ^Hl) depends on factorization 
in the ring of Gaussian integers. By making direct use of this idea we not only find 
independent generators for S'OC'(Z^) but also have a method that readily generalizes 
to other lattices and modules. 

To this end, we consider the lattice as the ring Z[i] of Gaussian integers, i.e., 
with i = y/ — 1, 

Z[i] = {a + ib I a,6 G Z} (11) 
together with the (number theoretic) norm 

noTm{a + ib) = {a + ib){a — ib) = \a + ib\'^ . (12) 

The ring Z[i] consists of all algebraic integers in Q{i) = {a + ib \ a,b E Q}, which 
is both a quadratic and a cyclotomic field. The coincidence rotation problem is 
then equivalent to finding all numbers of norm 1 in Q{i) because rotation through 
an angle means multiplication with the corresponding complex number on the unit 
circle and a coincidence can only happen if this complex number is in Q{i). 

Any such number can uniquely be written (up to units) as the quotient of two 
Gaussian integers, 

e- = ^ = ^ (13) 
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with coprime Gaussian integers a, (3 of identical norm, norm(a+i6) = norm(c+2(i) = 
say. Now, we can profit from unique factorization in Z[z] because any integer 
a G Z[i] divides its norm: 

a I norm(a) . (14) 

Let ^ = p'^ ■ ■ ■p^'^ be the (unique) factorization of I of (|13|) into "ordinary" primes of 
Z, called rational primes from now on. If any rational prime pj stayed prime in Z[z] 
(i.e., did not split into two Gaussian integers), which happens if pj = 3 (4), it would 
appear both in the numerator and the denominator of (|T3|l which is inconsistent 
with coprimality. Thus such a rational prime cannot divide £. 

A similar argument applies to the prime 2 which, although it splits as 2 = 2(1— i)^, 
also has only one Gaussian prime factor up to units. The remaining primes are = 1 
(mod 4) and split as p = Ujpjp into two Gaussian integers. One of them appears in 
the numerator of ()13|). the other in the denominator, ii p\ C.. Of course, the actual 
choice of Up is only unique up to units and up to taking the complex conjugate 
which reflects the point symmetry of the square lattice! One convenient choice for 
uniqueness (which we will now take) is a rotation angle in the interval (0, vr/4). 

This, in fact, solves the above problem constructively: any CSL rotation can be 
written in the form 

'■^=- n i^)" (15) 

V3p=l (4) ^ 

where G Z, e is a unit in and V denotes the set of rational primes. Since the 
group of units in Z[z] is nothing but C4, we find 

SOC{Z^) ~ C4 X Z(^«) (16) 

and the generators are i (for C4) and ujp/oJp for rational primes p = 1 (4). By 
we mean, as usual, the infinite Abelian group that consists of all finite integer 
linear combinations of the (countably many) generators. The coincidence index m 
is obviously 1 for the units in C4 and p = norm(ci;p) for the other generators because 
this counts the number of residue classes of the CSL in 7?. If the CSL rotation R 
is factorized as in Eq. |T3|). we thus find 

S(i?) = Yl pl'^^'l . (17) 

V5p=l (4) 

This solves the rotation part in principle, one can now work along the primes 
p = 1 (4) to write down the generators explicitly, e.g., 

4 + 3i 12 + 5i 15 + 8i 21 + 20i 35 + 12i 40 + 9i 

5 ' 13 ' 17 ' 29 ' 37 ' 41 ' ^ ^ ' 

where the number on the unit circle is shown in a form with denominator p and 
rotation angle in (0,71/4). All other CSL rotations are obtained by combinations, 
and one can regain the formula of Theorem Q for the number of them with index m. 
Since d*{m) is a multiplicative function (i.e., d*{mim2) = d*{mi)d*{m2) for coprime 
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mi, 7712) and d*{p'^) = 2 for a prime power p"^ (r > 1), we obtain for /(m) = d*{m) 
the Dirichlet series generating function jTH] 



oo 



d*{m) 



T-r / 2 2 



) 



m=l 




(18) 



and the Dirichlet series generating function for f{m) is 4$(s). 

Finally, the full group of CSL isometries, OC(Z^), is the semidirect product of the 
rotation part S'OC(Z^) (normal subgroup) with the group Z2 generated by complex 
conjugation (= reflection in the x-axis): 



Here conjugation of a rotation through an angle (p by complex conjugation results 
in the inverse rotation through —(f. Let us give a brief justification of Eq. ()19j) . 
Since 0(2) = 5*0(2) x^, Z2 (semidirect product) with the Z2 of Eq. (fT^ . any planar 
isometry T with det(T) = — 1 can uniquely be written as the product 



of a rotation through (f with T^., the reflection in the x-axis. But leaves 
invariant, so T is a coincidence isometry if and only if R{if) is a coincidence rotation. 

The calculation of coincidence indices is also simple in this case. The coincidence 
index for the reflection is 1. For an arbitrary element of OO(Z^), we either meet 
a rotation (where we know the result already) or use the factorization (|2(jp again. 
Then, the coincidence index is identical with that of its rotation part, so Eq. (|2(jp 
is all that is needed. This solves the coincidence problem for the square lattice 
completely and we have 

Theorem 2 The group of coincidence isometries of the square lattice 1? is 



This group is fully characterized by Eqs. (fT3j) . (f7^ and (f7^. and the coincidence 
index of an element is given by Eqs. ( [77| ) and ff73| . 

3 More generality: the unique factorization case 

As briefly explained in the introduction, the corresponding programme for a locally 
finite tiling T with A^-fold symmetry (or rather for its set of vertex sites) consists 
of two steps, the first being the solution of the coincidence problem for the limit 
translation module ^A.{T) of T (see [201 fo'^ details about this concept). For the 



00(Z2) = SOC{Z^) x,Z2 . 



(19) 



T = R{^) ■ T,. 



(20) 



OO(Z^) ~ 0(2, Q) ~ (O4 X Z(^«)) X, Z2 . 



(21) 
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moment, we consider only tilings with the property that the set of vertex sites of T 
is a subset of Ai{T), a condition we shall come back to in Section 5. Furthermore, 
we assume Ai{T) to be what is termed an 'W-lattice" in [2Tj but which we shall 
call an "n-module" (where N = 2n for n odd and N = n otherwise) in line with the 
mathematical practice of reserving the word "lattice" for discrete subgroups. The 
principal n-module (the "standard A^-lattice" of [21]) is the additive subgroup of M^, 
generated by the vectors of the regular n-star, 

(cos(27rA;/n),sin(27rfc/n)) , A; = 0, . . . , n - 1 . (22) 

The other modules are the non-trivial subgroups of the principal module that are 
invariant under rotation about the origin through 2tt /n. Modules that differ only in 
scale and orientation are regarded as equivalent. 

Because all modules are invariant under rotation through vr (since if x is in the 
module then so is —x), an n-module with n odd is invariant not only under rotation 
through 27r/n but also through vr/n. So n- modules and A^- modules are the same. 
In view of this we shall assume throughout that n is either odd or divisible by 4, 
though this necessitates bearing in mind that for odd n an n-module has 2n-fold 
symmetry. The opposite convention is used in [21], but the one used here is more 
convenient for expressing results about cyclotomic fields that we shall need later 
because it gives n the parity of the discriminant of the corresponding field. 

The first stage of our analysis, occupying all but Section El is to investigate 
coincidence rotations for modules and their associated coincidence site modules, 
which we designate CSMs. 

Symmetric modules and cyclotomic fields 

Viewed as complex numbers, the vectors ()22|) are where ^ is a primitive nth 
root of 1, and the modules are subsets of the cyclotomic field K = Q(^). The 
principal module is precisely the ring of integers Ok of K, since it is known that 
{1, ^, . . . , ^''^("')~^)} is a basis for the integers of K, where 0(n), the Euler totient 
function of n, is the degree of K over Q, cf. Chapter 9 of |22^. The other modules 
are the ideals of Ok (to be defined later), modules being equivalent precisely when 
they belong to the same ideal class (defined in Section E)). 

In this section, at the expense of discussing only 29 modules (see j^Hl HI])) we 
restrict attention to values of n for which all n- modules are equivalent. Because 
of the connection with algebraic number theory we call this the "class number 1" 
case and use the designation "CNl" to indicate results that are special to this case. 
(The reason behind this terminology is explained in Section 6. Briefly, it is the case 
when the nth cyclotomic field has class number 1.) The class number 1 assumption 
simplifies the treatment in two ways: 

1) it is enough to solve the coincidence problem for the principal module Ok 
only, since all others are equivalent to it; and 
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2) in the class number 1 case each integer in Ok has a factorization into irre- 
ducible integers that is unique apart from multiplying the factors by units. 
(Because of the unique factorization these irreducible integers can safely be 
called primes in the class number 1 case.) 

Though a convenience, the restriction to class number 1 is by no means essential: 
with only minor modifications our method applies to any 2D module, as outlined in 
Section El 

As in the previous section, a coincidence rotation that takes /3 to a, say, {a, (3 G 
Ok) can be represented by the point •y = a/P on the unit circle. So the CSM 
problem amounts to finding the structure of the set of numbers j in K with 

l7l = 1 (23) 

(a subgroup of the multiplicative group of K). 

The CSM associated with 7 is Ok H 'jOk = num(7)(9ii', where num(7), the 
numerator of 7, is given by 

num(7) = gcd(z/ e Ok \ «^/7 G Ok) , (24) 

and is unique up to multiplication by a unit. In particular, num(7) | a. The index 
of this module in the original module Ok is norm(num(7)), the absolute norm of 
num(7), (^j 4.4 and Cor. 2.96). (Since units have norm 1 this is independent of the 
particular numerator chosen. All conjugates of the field K are complex, so norms 
of numbers in K are products of pairs of complex conjugates and hence positive.) 

Eq. fl23|) can be reformulated as an algebraic condition with the aid of the maximal 
real subfield L oi K: 

^:=Q(e + r')=Q(cos— ). (25) 

n 

It is known that when K has unique factorization L does too (see p. 231 of |2Z1)- 
As an extension of L, K has degree 2 and the set of conjugates over L of a number 
7 G -ft' is just the complex conjugate pair {7,7}. Consequently, the relative norm 
of 7 over L, noTm.K/L{l), is given by 

nOTTClK/Lil) = ■ (26) 

In this notation, the absolute norm of 7 is norm(7) = normi^/Q(7) and we have the 
relation 

normi^/Q(7) = normL/Q(norm/^/L(7)) . (27) 

Relative norms of integers in K are integers in L and norms of units are units. As 
in the previous section (where L = Q), a \ noTmK/L{<^) = oia for every integer a of 
so the only possible prime factors of a in are those that divide norm(a). 



9 



Cyclotomic numbers on the unit circle 

When a planar module M. intersects a rotated or reflected copy of itself in a submod- 
ule of finite index, the isometry (rotation or reflection) is again called a coincidence 
isometry. The set of coincidence isometries of Ai is denoted by 0C{A4). It is again 
a group, with SOC{Ai) being its subgroup of rotations. (These concepts can be put 
in a much more general setting. Some slight extensions of them are already required 
for the examples in Appendix A, for example.) 

In view of Eq. (j26|) and our representation of SOC{Ok) as the elements of K on 
the unit circle, we have 

SOC{Ok) - {leK \ normK/Lil) = 1} • (28) 

To analyze the right hand side further we need some facts about the arithmetic 
of K and L. First, the units e of K with \e\ = 1 are precisely the powers of C,, 
though in general there are also infinitely many units not on the unit circle. (This 
follows, e.g., from j2Hl! Lemma 1.6, and the last sentence of the remark following 
it.) Second, if a prime g of L has two non-associated prime factors in K (i.e., their 
ratio is not a unit) then they can be taken as complex conjugates, u and uJ. This 
is because u \ q implies uo \ g, so, if u and ZJ are not associates, uZJ is an integer in 
Ol dividing g, hence is an associate of g. (Here Ol is the ring of integers of L, of 
course.) Conversely, if g is divisible by just the prime u in K and no other (up to 
units), then, as uJ also divides g, u/U must be a unit. Thus a prime u G Ok divides 
a prime g G Ol with distinct factors if and only if to is not an associate of cu. By 
Eq. (j2ni), noimK/Li^) = noruiK / l{oJ) . 

Now suppose that •y E K satisfies normK/Lil) = 1 and write 7 = where a, 
j3 are integers of Ok with no common factor. Then 

normx/L(a) = norm^/i(/5) = z/ G Ol (29) 

and every prime factor of v must factorize into two non-associated primes of K, one 
of which divides a only and the other (5 only. Since any such pair can be chosen to 
be complex conjugates, 7 can be written as 

-^n(l) 

with e a unit of K and the n^'s in Z. Taking absolute values in (jHUj) shows that 
\e\ = 1, whence e is a root of unity. Different values of the ra^'s give 7's with different 
prime factorizations, which are therefore not associates, and different roots of unity 
e give different 7's within each set of associates. So in this more general situation 
we again have explicit presentations of SOC{Ok) and OC{Ok) almost identical to 
those for S0C{7?) and 0C{7?) in the previous section. These are, for SOC, 

SOC{Ok) ^ (0 X ( = ) ^ C;vxZ(^«), (31) 
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(30) 



where Q is the set of complex conjugate pairs of non-associated primes in K (and 
N = lcm(n, 2) as usual) and, for OC, 

OC{OK) = SOC{OK)y<,{-r), (32) 

where ~ is complex conjugation and its action on SOC{Ok) is clear. 

The coincidence index of the typical rotation (jHUI) of SOC{Ok) is the absolute 
norm of its numerator: 

H (norm^/Q(cUfc))l"'=l . (33) 

k 

For each prime pair {uJ,uJ} G Q, the common value norm(u;) = norm(aJ) is 
a rational prime power p'^. We call these prime powers the basic indices of Ok 
and the primes p themselves the complex splitting primes for K (because, in their 
factorization over K, they contain at least one complex conjugate pair of distinct 
primes). Then (jSH), and show that: 

Proposition 1 (CNl) An integer m E N is a coincidence index if and only if it is 
a product of basic indices. 

To find out what basic indices there are and count how many members of the 
group SOC{Ok) have given index we need to determine how each rational prime p 
factorizes in the fields L and K. It will turn out that the basic indices are powers of 
distinct primes and that whether a power of p is a basic index (and what this power 
is) depends only on the residue class of p mod n. 

Factorization of primes in algebraic number fields 

Before considering K and L specifically we describe how primes factorize in a general 
algebraic number field extension F{a) D F of degree D. (We shall use the standard 
notation F{a)/F to denote such an extension. A detailed account of the material 
in this section can be found in Chapter 2 of 123)-) 

Let f{x) = be the minimal equation satisfied by a with coefficients in F (so 
the degree of f{x) is D) and let O and O' be the rings of integers of F and F{a). 

An ideal of (9 is a subset a of O (non-empty and 7^ {0}) such that a + (3 E a, 
Va, P E a, and Xa E a, \/X E 0,a E a. We use the notation (a, /3, . . .)o to denote 
the smallest ideal containing a,(3, . . . (where these are numbers in O). Ideals have 
a natural multiplication, defined by ab = [ajS \ a E a, (3 E h)^ and O itself is the 
multiplicative identity. There is an infinite set of prime ideals in O and every ideal 
can be uniquely factorized into prime ideals. 

Every ideal a in C extends to an ideal (a,)o' in C, but O' also has other ideals 
not of this form. For an ideal a' of O' the relative norm, noim.p(^a)/F{0'')-, is defined 
as {jioTniF{a)/F{.Oi) I a E a')o — an ideal of O. Norms are completely multiplicative 
(i.e., norm(a6) = norm(a)norm(6)). Let p be a prime ideal in O. Then {p)oi 
factorizes into prime ideals in O' as 

(p)o'=Pr...p^% (34) 
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and for each k = 1, . . . ,g 

norm^(„)/p (Pfc) = P"^' , (35) 

where dk is the residue class degree of p^. Taking norms of both sides of Eq. 
shows that 

diCi H h dgCg = D, (36) 

For the special case of normal field extensions (i.e., extensions where F{a) contains 
not only a itself but also all other roots of f{x) = 0) we have Ci = ■ ■ ■ = Cg and 
di = ■ ■ ■ = dg. So in this case 

(p)o' = (Pi...p,)% (37) 

where each has the same degree d, d ■ e ■ g = D and without ambiguity we can 
define degp(^^^{p) := d and e^(.^-)(p) := e. Also Gp^^-^i^p) > 1 only for the finitely 
many primes p that divide the discriminant of the extension F{a)/F. (Such primes 
are called ramified.) 

Another special case is field extensions where a can be chosen so that O' = 0[a] 
{O' has a simple integral basis over O). In this case the factorization of a prime p 
of O into prime ideals of O' mimics the factorization of f{x) into irreducible factors 
over the finite residue class field O/p. So if p factorizes as in Eq. (jMjl then 

fix)^Mxr-..fgixr^ (modp) (38) 

where each fk is irreducible of degree dk and distinct /a,.'s correspond to distinct 
primes p^- This provides a simple way of calculating the degrees and multiplicities 
of the prime factors of p. 

The three extensions we have to deal with — K/Q, L/Q and K/L — are all 
normal and have simple integral bases, so all the above results apply to them. Also 
relative degrees are multiplicative: if p is a rational prime having a prime factor g 
in Ol which in turn has a prime factor u in Ok then 

degi^/Q(a;) = degj^/L{uj) -degL/Qig). (39) 

In particular, deg^/^(cij) is the same for all prime factors u in Ok of the same 
rational prime p. 

The primes of Ok in the non-associated pairs {uj,uj} are precisely the unramified 
primes of relative degree 1 over O^. In view of (I39p and the normality of K and L, 
Q is the set of all pairs of distinct prime factors {uj, uJ} in Ok that divide rational 
primes p with 

degL(p) = deg^(p) (= d, say) (40) 
and, for any such u, the absolute norm (cf. ^I7\i above) is 

norm(a;) = p'^. (41) 

We have: 

Proposition 2 (CNl) The complex splitting primes for K are the rational primes 
that satisfy ^^'^ ^he basic indices of Ok are the powers p'^ of these primes. 
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How to calculate the CSL group and its coincidence indices 

Getting more explicit information about SOC{Ok) and its coincidence indices comes 
down to finding degi{p) and deg^(p) for rational primes p. The following facts are 
sufficient to do this; we state them here and justify them in Appendix B. 

To reiterate our notation: K = Q(^) and L = Q{C, + C,~^), where ^ is a primitive 
nth root of 1, and p is any rational prime. 

Fact 1 If p\n then deg;^(p) is the smallest (i G N such that n divides p'^ — 1. 

Fact 2 If p \ n then deg^(p) is the smallest d G N such that n divides at least one 
of p'^ + 1 or p'^ — 1. 

Fact 3 (a) If n = p^ , for some r, then p is not a complex splitting prime and 
degxip) = 1- 

(b) More generally, if n = p^ui with p \ ui then p is a complex splitting 
prime in K if and only if it is a complex splitting prime in Ki ( the cyclotomic field 
of Tilth roots of unity). Moreover, deg^(p) = deg^^(p). 

Although these facts alone clearly enable us to identify the complex splitting 
primes and calculate their degrees and multiplicities, it is nevertheless worth listing 
some general consequences of them. 

Remark 1 These facts show that whether a prime p is a complex splitting prime 
of n and what its degree d is depend only on the residue class of p mod n. 
Remark 2 Since, for a prime uj in Ok dividing p, norm(c(j) = where d = 
degi^{p), Fact m has the well-known consequence that norm(co') = 1 (mod n) for 
every prime u in Ok with u \ n. In particular, every coincidence index m with 
gcd(m, n) = 1 satisfies m = 1 (mod n). 

Remark 3 When p \ n is not a complex splitting prime, deg^(j9) = 2 degi^{p), so 
degx(p) is even. Fact s[TJ|21 and 3(a) show that, conversely, if n is an odd prime power 
then no prime p with d = deg^(p) even is a complex splitting prime. This is because 
if n I = (p'^/2-l)(p°'/2_^l) but n f p"'/^ _ X then, since gcd(p'^/2-l,p'^/2 + l) = 2 

(or 1 if p = 2), n I + 1, so deg^(p) = d/2. 

So, for n an odd prime power, p is a complex splitting prime if and only if degj^(p) 
is odd and p\ n, and it is unnecessary to compute degrees over L in this case. 
Remark 4 By Fact the unramified primes with deg^(p) = 1 (i.e., the primes 
that split completely in K) are precisely those = 1 (mod n). So these primes are 
always complex splitting primes. 

Remark 5 Facts Hand 121 show that, for primes p = —1 (mod n), degj^(p) = 2 and 
degi{p) = 1, so these primes are never complex splitting primes. Consequently, for 
every n, the proportion of integers that are coincidence indices is 0. 

In the next section we apply these facts and remarks to calculate coincidence 
indices of specific modules. 
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^ (42) 



The number of coincidences with given index 

Let fivn) = N ■ f{m) be the number of elements of SOC{Ok) with index m. The 
computational convenience of representing f{m) this way arises from the fact that 
f{m) is more fundamental: it is a multiplicative function of m and, as for the square 
lattice, it counts the CSMs with index m, since the rotation group of each module is 
the group of roots of unity in K and has order N. In the general case, f{m) cannot 
be described as simply as in Eq. ©, but its Dirichlet series generating function does 
have a very simple expression in terms of the (^-functions of the fields K and L. Also, 
for quite sizeable individual values of the index, the number of coincidence isometrics 
can be calculated from (jH^ and knowledge of the identity, degrees and exponents 
of the complex splitting primes (or, equivalently, from the generating function). 

From the decomposition of SOC{Ok) and the function (jH^ it can be seen 
that /(m) is multiplicative, i.e., gcd(mi,m2) = 1 implies f{mim2) = f (mi) f {1712) ■ 
This makes its Dirichlet series 

fjjn) 

m=l 

a convenient tool for studying /: it can be expressed [IHj as an "Euler product" 

n(l^). (43) 

with one Euler factor for each prime p, and the individual Euler factors are straight- 
forward to compute. (The series we obtain will all be absolutely convergent in the 
right half-plane Re(s) > 1 and extendable to meromorphic functions on the whole 
plane. For using the series formally to calculate individual values of / these an- 
alytic properties are irrelevant, but they play an essential role in calculating the 
asymptotic average value of /.) 

Suppose the rational prime p is divisible by the pairs {cji, cJi}, . . . , {ujg/2, ojg/2} of 
non-associated primes in K and that each ujj has noivuK/Qil^j) = p'^- Then f{p^) is 
the coefficient of p~^^ in 

11 11 y/^ 

~^ p2ds ~^ pds ~^ ~^ p(ig ~^ p2ds ~^ ) ' (^^) 

the product of g/2 two-way infinite sums (one for each pair {ujj,ujj}) each having 
one term for each value of the corresponding Uk in ()33|1 . (The symmetry of the sums 
arises from the fact that the index depends only on \nk\, of course.) On summing 
the series this becomes 

Since f{m) is multiplicative, for a general m it is the coefficient of in 

n (1^)"^ 
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where C is the set of complex sphtting primes for K, and the values of d and g are 
those appropriate to each individual prime p. This, in turn, is the coefficient of m~* 
in the infinite product 



n 

pec 



l+p- 
1 - p-'^" 



ds \ 9 12 



(47) 



To express this more simply we introduce the Dedekind ^-functions of number 
fields j22t I2H1 • The (^-function of a general algebraic number field F is the Dirichlet 
series generating function for the number of ideals a oiO with norm(a) = m, hence 
is given by 

c^(^) = E — W = n — W) ' (48) 

norm(a)'' p V norm(p)*y 

where a runs through all ideals of F and p through all prime ideals. When F is 
normal we can collect together prime ideals p dividing the same rational prime p to 
put the product on the right in the form 



J. J. pds 



(49) 



where, for each rational prime p, d = degp/Q{p) and g is the number of prime ideals 
of F dividing it. 

A particular case of this is 



oo ^ 



m=l 



1 

pS 



-1 



(50) 



which is the Riemann ^-function ({s) itself. 

The following table compares the Euler factors of Ck{s) and Cl(2s) for each 
rational prime p, there being three cases to consider. (It follows from Lemma 3 of 
and Prop. 2.15(b) of ^\ that the third case, exip) 7^ ^l{p), occurs for at most 
one prime p: the prime, if any, a power of which is equal to n.) 



p 


Field 


Degree 


Distinct 
prime factors 


Euler factor 


Complex splitting 


K 
L 


d 
d 


9 
9/2 


(1 -p-'^')-9 
{I _ p-2dsyg/2 


Not complex splitting 


K 


d 


9 


(1 -p-'^')-9 


and exip) = eiip) 


L 


d/2 


9 


(1 -p-'i')-9 


exip) + Clip) 


K 
L 


1 
1 


1 
1 


{l-p~')-^ 
(1 _p-2s)-l 
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On taking the quotients of the Euler factors arising from K and L and comparing 
with Eq. (jlTj), we see that 

the star indicating that the second factor on the right is present only if n is a power 
of a prime p. So 

( {1+P^-')-^Ck{s)/Cl{'2s), if n is a power of a prime p, 
CA'l =^ '^k[s) = \ (52) 
\CK{s)/(Li2s), if not. 

We summarize this in 

Theorem 3 Let n he one of the 29 numbers for which the cyclotomic field K of 
nth roots of unity has class number 1. Then the group of coincidence rotations of 
an n-fold symmetric module is the direct product of its finite rotation symmetry 
group Cn o,nd countably many infinite cyclic groups, as in . and the full group 
of coincidence isometrics is the extension of this by a reflection symmetry. The 
coincidence index of such an isometry is given by 1^3^ and (jj^ . The Dirichlet 
series generating function for 

f{m) = {number of CSMs of index m} (53) 
= ^ {number of coincidence rotations of index m} (54) 



IS 



given by (t5gjl . 



A principal use of a Dirichlet series is to find asymptotic formulae for sum func- 
tions of its coefficients by means of residue calculus. In the present instance this 
technique shows, for example, that 



Number of CSMs of index <X = ^ f{ 



m<X 

~ X ■ {residue of ^k{s) at s = 1}. (55) 

In view of Eq. (jHH), this residue can be computed from known formula for the 
residues of (^-functions at 1 and values of (^-functions at 2. The value of the residue 
can be regarded as the "average number of CSMs" with a given arbitrarily chosen 
positive integer as index. 



4 Examples: iV = 6, 4, 10, 14, 8, and 12 

After the general derivation of the previous section, let us present some examples 
explicitly. We select those relevant to known crystals and quasicrystals. For each 
example we list: 
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(a) The fields K and L and the degree [K : Q] of ii" over Q. 

(b) A table giving, for each residue class mod n containing primes p, deg^(j9) (and, 
if necessary, deg^(p) too). In the bottom line of the table (where deg^(p) is 
given) the degrees of complex splitting primes are underlined. With each table 
is a comment describing which facts and remarks from the previous section 
were used to compute it. 

(c) A list of the types of basic indices, using the notation that p^^^ represents the 
6th powers of all primes congruent to a mod n. 

(d) The Dirichlet series generating function of f{m), given as a ratio of (^-functions, 
as an Euler product and expanded explicitly as far as the 12th nonzero term. 
(The same notation as in (c) is used for the primes in the Euler product.) 



(e) An explicit formula for /(m) in the style of Eq. 
denotes the largest exponent e for which | m. 



In these formulae Cr. 



(f) The average value of f{m), as defined above. 

The smallest coincidence indices can be read off as the denominators (with s = 1) 
of the Dirichlet series, with the corresponding values of f{m) as the numerators. All 
values of f{m) for m > 1 are even, reflecting the geometrical fact that the reverse 
of a coincidence rotation is also a coincidence rotation. 



n = 3, the triangular (or hexagonal) lattice 

ir = Q(v^), L = Q, [K:Q]=2. 



p (mod 3) 


1 2 3 


deg^(p) 


1 2 1 



Computed using 

Facts [Hand 3(a) and Remark 3. 



Basic indices: p^^^-^ 
Dirichlet series: 

-1 



1 



Of(£) 
C(2s) 



Hi 



1) 



(1) 



ii-7Fi-Y3Ti-Yp-i-3pi-37Ti-4p-i-49Ti-gpi-g7ji-73?i-7p--|- 



Number of CSLs with index m: 

/(m) = 

Average number of CSLs: 



J^2, if m is a product of basic indices, 
0, otherwise. 



^ ^ 0.276 
27r 
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n = 4, the square lattice 

K = Q{^,L = Q, [K:Q] = 2. 



p (mod 4) 


12 3 




1 1 1 
1 1 2 



Computed from 
Facts d Eland 3(a). 



Basic indices: p^-^^^ 
Dirichlet series: 



Ck{s) 



Hi 



1) 



C{2s) --i-p^^ 

+ d 

Number of CSLs with index m 
Average number of CSLs: 



1 + — -I — - — I — - — I — - — I — - — I — - — I — - — I — - — I — - — I — - — I — - — I- 



J^2, if m is a product of basic indices, 

p\m 

0, otlierwise. 



1 



0.318 



n = 5, the 10- fold module 

K = Q{e^-'^/^), L = Q{^/5), [K : 



p (mod 5) 


1 2 3 4 5 


deg;^(p) 


1 4 4 2 1 



Computed using 

Facts HI and 3(a) and Remark 3. 



Basic indices: p^^^^^ 
Dirichlet series: 



1 



Hi 



1) 



p 



(1). 
4 



Cl{2s 

^ 31" 

Number of CSMs with index m: 

/(m 

Average number of CSMs: 



J- -r l]_s -r 3;^s T T g]_s -r 7;^^ t ^q^s t ^21" ^ 131" ^ I5l" ^ 181" ^ 191" ^ 



JJ^ 4ep, if m is a product of basic indices, 
0, otherwise. 



5 logr 



0.244 
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n = 7, the 14-fold module 



K 



L = Q(cos(27r/7)), [K : 



p (mod 7) 


1 2 3 4 5 6 7 




1 3 6 3 6 2 1 



6. 



Computed using 

Facts Hand 3(a) and Remark 3. 



Basic indices: p^^^^^, p^2)^ pfi) 
Dirichlet series: 



1 + 



1 Ck{s) 



a(2s) 



6 



1 _|_A4__6__|__6__|__2__|__6__| 6 I 6 , 6 , 6 , 12 6 , 

8= 29= 43» 64^= 71" 113'' 127= 197" 211» 232" 239** 



Number of CSMs with index m: 

Y\ (4ep + 2) Y\ 2, if m is a product of basic indices, 



/(m) 



p|m 
p=l (7) 



p\m 
P^l (7) 



otherwise. 



Average number of CSMs: 

where R (the regulator of K) is given by 



0.235, 



J = log2(2 cos y) - log(2 cos log(2 cos y) ^ 0.525. 

n = 8, the 8-fold module 

K = Q(e-/^), L = Q{V2), [K : Q] = 4. 



p (mod 8) 


1 2 3 5 7 


degi(p) 
deg^(p) 


112 2 1 
112 2 2 



Computed from 
Facts mil and 3(a) 



Basic indices: p^^^, p^g-j, p^^^ 
Dirichlet series: 



Cl(2s) 



Hi 



■L -r gs -r ;^73 -r 25s ^ 41s ^ 73s ^ g^s T ggs -r gys T ^^^s T 



121' 



+ 



137' 



+ • 
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Number of CSMs with index m: 

Y\. n ^' if is a product of basic indices, 



/(m) 



p=l (8) p^l (8) 



0, otherwise. 
Average number of CSMs: 

2x/21og(l + y2) 



0.253 



n = 12, the 12-fold module 



K 



, [ir : Q] = 4. 



p (mod 12) 


1 2 3 5 7 11 




1112 2 1 
1 2 2 2 2 2 



Computed from Facts 1 and 2 
and the cases n = 3 and n = 4 
using Fact 3(b). 



Basic indices: p^-^^ pj^y p^j^ 
Dirichlet series: 



Cl{2s) ^^[l-pl,]) (1-P(5f)(l-P(,f) 

-L -r 13s 25" ^ 37" ^ 49" ^ 61" ^ 73" ^ 97" ^ 109" ^ 157" ^ 169" ^ 181" ^ 

Number of CSMs with index m: 

Y\. n ^' if is a product of basic indices, 



/(m) 



p\m p\m 
p=l (12) p^l (12) 



0, otherwise. 
Average number of CSMs: 

y3Iog(2 + v^) 



0.231 



5 Application to 2D quasicrystals 

The reader might like to see at least one or two examples where we apply the above 
results to planar quasicrystals. For simplicity, we consider the eightfold symmetric 
Ammann-Beenker tiling and the decagonal Tiibingen triangle tiling [22] here, while 
the slightly more complicated rhombic Penrose tiling is discussed in Appendix A. 
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Figure 1: Central patch of the exactly eightfold symmetric Ammann-Beenker tiling. 
The Ammann-Beenker tiling 

Consider the eightfold symmetric Ammann-Beenker tiling of Fig. 1 and, in particu- 
lar, the coincidence problem of its vertex points for rotations around the symmetry 
centre. The underlying module is the standard eightfold module of rank 4, usually 
obtained as projection of the hypercubic lattice to a suitably chosen 2D plane. 
This plane, and its perpendicular complement, are eigenspaces of an eightfold rota- 
tion. 

The set of vertex sites of this tiling is just the subset of module points whose 
corresponding points in perpendicularly project into a certain regular octagonal 
window. It is clear then that a coincidence of vertex sites implies one in the module, 
but also the converse is true due to the way the tiling sites are distributed over the 
module. 

A coincidence rotation can be hfted to 4-space whence it also affects the window. 
In fact, a coincidence point must have perpendicular projections both in the original 
and in the rotated window! But this results in a slight modification of the fraction 
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Figure 2: Intersection of two acceptance domains that are rotated against each other. 



of coinciding points which has to be corrected by an acceptance factor A. This 
is nothing but the area ratio of the intersection of the rotated windows with the 
original window, see Fig. 2. For a coincidence rotation through ip, it turns out to be 

A = 1 - (1 - i=) sin(V;) sin(^ - i^) , (56) 

where G [0,7r/4) via 

ijj = ilj- — ■ - 57 

L TT J 4 

and ip, the rotation angle in perpendicular or internal space, is related to the angle 
Lf = 2arctan(a + b^/2) through an algebraic conjugation: 

ip = 2 arctan(a - bV2) . (58) 

The acceptance factor (j^UI) is 1 for symmetry rotations and smaller otherwise, 
the minimum value being Amin — 0.957 at tt/8. The set of coinciding points almost 
looks like an Ammann-Beenker pattern again, but some points are missing: the 
quantity 1 — A is the frequency of such failures which were observed in j^j. With a 
more complicated window, star-shaped say, the acceptance factor would also become 
more complicated: with some choices of window it can even be zero for certain angles. 
But we will not go into further details here. 

The Tiibingen triangle tiling 

Let us now consider the coincidence problem for the vertices of the decagonal trian- 
gular tiling of Fig. 3. All vertex sites belong to the standard tenfold module which 
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can be obtained by projection of the root lattice A4 to a suitably chosen plane 
For simplicity, we consider the cartwheel tiling (which is singular) because it has 
full Dio symmetry in the sense that a Diq operation produces mismatches of density 
zero in the plane (along worms). We thus have coincidence fraction 1 in this case. 
Also, all other coincidences of the tenfold module are realized. As in the previous 
example, one has to correct the coincidence fraction, this time by rotating a decagon 
(the window of the vertex sites) and intersecting it with the original one. Let us 
give the correction formula in slightly more generality. If the window were a regular 
77,-gon, the analogue of Eq. would read 

/sin(Q;/2)\^ , , , 

A=l - { — Vt^ sinV' sin a - tp) , 59 
\ sm(a) / 

where a = 2tx /n, ip = — ~, and ip is related to ip via an algebraic conjugation. 

In the present instance, this relation is that tan(y9/2) can be expressed in the 

form 2 

tan ^ = (a + hr) sin — (a, & G Q) (60) 
2 5 

and then 

tan — = (a + or ) sin — , (61) 
2 5 

where r' = — 1/r is the conjugate of r in Q(t). 



Twelvefold symmetric tilings 

As well as eight- and tenfold symmetries, twelvefold symmetry is of practical interest. 
Here the calculation of ip is very similar to the eightfold case: given the angle 

= 2 arctan(a + hV?>) (62) 
in tiling space one obtains the angle 

ip = 2 arctan(a — h\/?>) (63) 
in internal space, which can be used with Eq. ()59|) . 



6 Beyond unique factorization 

In Section 3, we restricted ourselves to the "class number 1" case, where there is 
essentially only one ra-module. We now show how our method can be adapted to 
other cases, too. The smallest value of n to which Section 3 does not apply is 23 
(A^ = 46), mentioned in [21]. Here, the cyclotomic field has class number 3, so there 
are 3 distinct modules with 46-fold symmetry. (The number of modules increases 
rapidly with n [211^1 •) 
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Figure 3: Central patch of the cartwheel version of the tenfold symmetric Tiibingen 
triangle tiling. 
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6.0.1 Ideals and ideal classes 

Let F be algebraic number field with ring of integers O. The set of ideals of O can 
be extended to form a group by admitting fractional ideals of the form 

ah ^ = {7 I 7/5 G a V/3 G 6} , (64) 

where a and h are ideals as defined in Section 3. (A fractional ideal need not be a 
subset of O.) The identity element of the group of fractional ideals is O. A principal 
ideal is a fractional ideal of the form 

(7)0 = {7« I « G O}, (65) 

generated by the single number 7. When unique factorization into irreducible inte- 
gers fails in F then some ideals must necessarily be non-principal. Two fractional 
ideals a and b are equivalent if 6 = 7a for some 7 G -F. The equivalence classes, 
called ideal classes, form a quotient group of the group of fractional ideals called the 
ideal class group, H = H{F), which turns out to be finite. Its order is called the 
class number, h{F). The identity element of H is the class of principal ideals. 

The ideal classes inherit complex conjugation from F: each ideal class C E H 
has a complex conjugate class C. Complex conjugation is an automorphism of H of 
order 2. 

6.0.2 Ideals as modules 

Our definition of n-modules makes them ideals in the ring of integers of the nth 
cyclotomic field (and with any broader definition an n-module would certainly be 
equivalent to one of these). Multiplication by a complex number 7 is equivalent to a 
combined rotation and scale change in the plane, so equivalent ideals certainly cor- 
respond to equivalent modules. Conversely, equivalent modules can be transformed 
into each other by multiplication by a complex number 7, and if both modules are 
subsets of an algebraic number field K then 7 is in and the corresponding ideals 
are equivalent. 

So the set of n-modules up to equivalence corresponds to the class group of the 
nth cyclotomic field. 

6.0.3 Coincidence rotations in the general case 

With class number > 1, ra-modules are no longer all equivalent. So, for compre- 
hensiveness, we need to consider not just OC{Ok) but also OC(c) for an arbitrary 
ideal c of Ok- 

There are two problems to be overcome in extending our method to the general 
n-module: 

(1) How to classify which of the products on the right of Eq. (|3(jp give rise to 
numbers 7 with I7I = 1 (when some cu^'s are non-principal ideals) and 
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(2) how to choose a representative of the reflection coset of OC, for modules not 
invariant under complex conjugation, and how to calculate coincidence indices 
of reflections from it. 

In this subsection we address the first of these. 

For the fractional ideal aa~^ to give rise to a number 7 G -ft' with I7I = 1 two 
conditions are necessary (and the conjunction of these conditions is also sufficient). 
They are 

(A) aa~^ is principal, and 

(B) for every 6 such that {6)0 = aa^, 66 = ee for some unit e of K. 

Condition (B) arises because 7 = e5 in Eq. (30) gives 66 = e^^e^^. Condition (A) 
is tantamount to saying that the ideals a and a are equivalent, in other words that 
a belongs to a class in Hi, the subgroup of H consisting of classes C with C = C. 
It is easily checked that Condition (B) also depends only on the class of a and is 
preserved under multiplication and inversion of classes. For Condition (B) to be 
applicable at all a must belong to a class in Hi. Consequently Condition (B) is 
equivalent to a belonging to a class in a certain subgroup H2 of Hi. 

When Condition (B) is satisfied the numbers 7 = C,66~^, where ^ runs through 
the roots of 1 in K, satisfy I7I = 1. In this case num(7) is the ideal a and can still 
be defined exactly as in Eq. (j211), provided that "gcd" is interpreted as meaning "the 
ideal generated by". Again c fl 7c = num(7)c and the coincidence index associated 
with the rotation 7 is norm(num(7)) (independent of the ideal c). In the general 
case, when the cu's may be non-principal ideals, a member of the product group on 
the right of Eq. (31) is a pair (root of unity, fractional ideal of the form aa~^) and 
the argument of Section 3 shows that elements of SOC{c) correspond precisely to 
those pairs with the class of a in H2. Such pairs form a subgroup of finite index in 
the full product group. One can choose a set of generators for this subgroup in much 
the same way as one chooses a basis for a lattice of finite index in a given lattice, 
and as in that case there is an infinite number of such bases and no canonical choice. 

Although SOC{c) has independent generators as a group, the set of coincidence 
indices in general no longer has independent generators as a semigroup. 

6.0.4 Coincidence reflections in the general case 

Our second problem was how to calculate indices of coincidence refiections for a 
module class in which no module is invariant under complex conjugation. Choose, 
for simplicity, a prime ideal p in the class (which is possible since every ideal class is 
known to contain infinitely many prime ideals). Then pHp = pp has index norm(p) 
in p. Every coincidence refiection of p has the form p = 7~ for some 7 G C with 
I7I = 1. Being a coincidence refiection on p, p{a) = f3 for some a,f3 G p. Hence 
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= j3 / a & K . The index of p is the index of p fl 7p in p which is 
norm(num(7))norm(p), if p /num(7), 

(66) 

norm(num(7))/norm(p), if p \ num(7). 

We note that there is a reflection of index 1 if and only if the class of p is in H2 
(when we can choose 7 to be a generator of the fractional ideal p/p) and that in that 
case dnni) agrees with our previous way of calculating the index. When the class of 
p is not in H2 the smallest reflection index is got by taking (7) = pa/pa, where a 
is the ideal of minimal norm such that the class of pa is in H2. Of course, OC{p) is 
OC{p) conjugated by reflection in the x-axis (corresponding isometries having the 
same index). We note that this is consistent with just replace p and 7 by their 
complex conjugates. 

Theorem 4 The group of coincidence rotations of a general n-fold symmetric mod- 
ule is the direct product of its finite rotation symmetry group Cn and countably many 
infinite cyclic groups which can be effectively computed and depend only on n. The 
index of any coincidence rotation so presented can be calculated explicitly. Any such 
module is equivalent to some prime ideal in the cyclotomic field of nth roots of unity, 
and in this form complex conjugation represents the coset of coincidence reflections 
whose indices can be computed from 1^) (they depend not only on n but on the 
individual module). Such a module need not have exact reflection symmetry. 

The following table list some statistics for the first few cyclotomic fields with 
h > 1. We follow Washington^^ in listing fields with their degree, 0(n), as the 
primary order and n as the secondary order. For each field we give n, A^, H, Hi, 
H2, the smallest rotation index and the smallest reflection index of the non-principal 
modules (for the principal module it is always 1). In brackets after each index we 
give the number of different rotations or reflections with that index. For all fields 
on our list H2 is the trivial subgroup consisting only of the identity element E of 
H. Also complex conjugation acts on the class group as multiplicative inversion for 
all these fields. 



Min. rotation Min. reflection index of 



n 


N 


Degree 


H 


Hi 


H2 


index 


non-principal modules 


23 


46 


22 


C3 


[E] 


[E] 


599 


(22) 


47 


(11) 


39 


78 


24 


C2 


C2 


{E} 


157 


(24) 


13 


(2) 


52 


52 


24 


C3 


[E] 


{E} 


313 


(24) 


13 


(1) 


56 


56 


24 


C2 


C2 


{E} 


64 


(2) 


8 


(2) 


72 


72 


24 


C3 


{E] 


{E} 


729 


(2) 


9 


(1) 


29 


58 


28 


CI 


CI 


{E} 


4931 


(28) 


59 


(4) 










[E] 


[E] 




(30) 


r 32 


(1): order 9 


31 


62 


30 




5953 


\ 125 


(5): order 3 



The two sets of figures in the last entry are due to the fact that non-principal 
modules with different orders in the class group of Q(e^'^*/'^^) have different minimum 
reflection indices. 
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6.0.5 Another example: = 46 

To illustrate the results of the previous subsection we treat in detail the case n = 23 
(with 46-fold symmetry). For this n, the class group H of K is H = {£,0,0"^}, 
where = E and C = C^. Hence Hi = {E} and therefore H2 = {E} too. 

The methods of Section 3 show that the complex splitting primes are precisely 
those that are quadratic residues mod 23 and for these deg(p) = 1 or 11 according 
to whether p = 1 mod 23 or not. The prime ideals p of Ok that divide a given 
rational prime p are either all principal or all non-principal (because the Galois 
group Gal{K/Q) permutes them transitively) and in the non-principal case fall into 
complex conjugate pairs of ideals, one from each of the classes C and C^. We 
partition the set of pairs Q into the sets Qi, Q2 as follows: 

ni = {{uji,uJi},{uj2,uJ2},...} (67) 

^2 = {{Pl,Pl},{P2,P2},- ■ ■} , (68) 

where the tUj's are numbers (corresponding to principal ideals) and where in Q2 we 
have chosen E C , p^ (z for each i. Finding all numbers of K on the unit 
circle is equivalent to finding all principal ideals with K/L-norm equal to Ol. (The 
numbers 7 are then the sets of associates of the generators of these ideals.) These 
ideals are precisely those of the form 

, i) n(l) 

with Uk divisible by 3 (since each Pf.p'^^ belongs to the class of order 3). This 
group of ideals has each uJi/uJi as an independent generator of the first factor, and a 
set of independent generators of the second factor can be chosen as follows: 

P1P2/P1P2, ■■■■ (70) 

Although this exhibits SOC{c) as having independent generators as a group, 
the set of coincidence indices no longer has independent generators as a semigroup. 
Instead of basic coincidence indices one has the prime powers 

p (p = l(mod23)) and (p = 2, 3, 4, 6, 8, 9, 12, 13, 16, 18 (mod 23)) 

which can be partitioned into two classes Pi and P2 (corresponding to Qi and ^2) 
according to whether or not the prime ideals dividing p are principal. As examples: 

599, 691, 829, 59^\ 101" G Pi 

and 

47, 139, 277, 461, 967, 3", 13", 29", 3l", 4l", 7l", 73" G P2. 

These examples were computed using the observation (derived from the last para- 
graph of Chapter 1 of [2H1) that p factorizes into principal primes if and only if it 
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factorizes into principal primes in Q(y/—23). A necessary and sufficient condition 
for this is the solubihty of the Diophantine equation Qx"^ + xy + y"^ = p. 
The general product of these numbers has the form 

m = • • ■p':^{pI],^T^+' ■ ■ ■ {pI'T^ X {Pi-factors}, (71) 

where Pi,...,pl^ are in P2 with = 1 (mod 23) and Pr+i,---,Ps ^ 1 

(mod 23). Now, for k = r + 1 . . . s define 

e.^!" '"I"" (72) 
[1 if not. 

Then m is a coincidence index if and only if 

ai H \- ar + er+i H h 7^ 1. (73) 

[The reason for this is that in choosing a principal ideal giving index m we can 
arrange that ^ is divisible by 3 in ()69|) by changing the sign of some n^s provided 
at least two n^s are not divisible by 3. Primes of degree 11 are divisible by only 
one pair of primes in K, but primes of degree 1 are divisible by 11 such pairs, so for 
these we can easily arrange that no rik is divisible by 3.] 

Consequently the ffist three rotation coincidence indices are 1, 599, 691, the 
smallest not composed entirely of primes = 1 (mod 23) is 2^^47 = 96256 and the 
smallest with no prime factors = 1 (mod 23) is 2^^3^^ = 362797056. 

The Dirichlet series generating function of /(m) can be found much as before, 
except that the contribution from non-principal ideals with K/L-norm equal to Ol 
must be omitted. This can be done using the three characters of the class group: we 
form three Dirichlet series (Hecke L-series), one for each character, by multiplying 
each norm in the series by the value of the character on its ideal. The required 
generating function is then the average of these three series. 

For the principal character (identically equal to 1) the corresponding Dirichlet 
series is exactly as in Eq. namely 

1 + — V'^. (74) 
23^ J a(2s) ^ ' 

For a non- principal character x the Euler factor for a prime p occurring in Pi is 
exactly as in (j44j) and (j45j) . For a prime p occurring in P2, however, the Euler factor 
is 

p2ds pds pds p2ds J I ]^ p—ds i ' ^ ' 

where 77^ = 1. Since this does not depend on which primitive cube root of unity 7] 
is, the Dirichlet series formed with the characters x x are the same and we have 
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/(m) 



m=l 
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where the first product is over the primes p in P2 (which are = 1 (mod 23)) and the 

second is over the 11th powers q in P2. 

In hne with ear her examples we give the first 12 nonzero terms: 

22 22 22 22 110 22 22 22 22 22 22 

599^ 691'' ^ 829^ ^ 1151'' 2209'' 2347^ ^ 2393^ ^ 3037'* ^ 3313" ^ 3359'* ^ 4463'' ^ " ' 

Note that this apphes to all three modules. 

For the principal module the reflection indices are the same as the rotation indices. 
For the two non-principal modules, however, the first three reflection indices are 47, 
139, 277. 



7 Concluding remarks 

Let us summarize our results. We have solved the coincidence problem for planar 
patterns with A^-fold symmetry by number theoretic methods. The first stage con- 
sisted of the analysis of lattices and modules in the plane where the coincidence 
indices are integers. 

For various cases of interest we have given the solution explicitly, in particular 
describing the set of possible coincidence indices and the number of coincidence 
isometrics with given index. The method is described in sufficient detail to allow 
other examples along these lines to be worked out. This is relatively easy for < 46, 
but the complication increases astronomically for larger as foreshadowed even in 
the example = 46, where the class number is only 3. 

The second stage was the explicit investigation of discrete structures associated 
with a given module. Here, in the non-periodic case, the calculation of the coin- 
cidence ratio requires a non-integral correction factor. We have demonstrated its 
calculation in several examples. 

Furthermore, the approach via algebraic number fields automatically yields sets 
of independent generators for the CSM group and therefore an explicit description 
of it. The group structure is interesting in itself because we deal here with infinite 
discrete groups that are countably generated and the structure of such groups is not 
at all obvious. 

An obvious next step is to extend the investigation to 3D examples. This is not 
only an interesting extension of the technique, but may have concrete realizations. 
There are two cases to consider: first the T-phases, i.e. quasicrystals which have a 
unique quasiperiodic plane and are periodic in the third direction. The CSMs for 
rotations around the unique axis are the ones treated in this paper. CSMs around 
other axes occur only when special relations hold between the lattice constants in 
the plane and perpendicular to the plane, a result familiar from the hexagonal case 
[21 . There are also near-coincidences with small misfits between the two grains, but 
it is beyond our scope to deal with these. The second case is the icosahedral one, 
the only remaining non-crystalline symmetry in 3D. Here we do not have such a 
powerful tool as the complex numbers and the structure of the CSM groups is more 
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complicated, even the rotation part being non-Abelian in general. Some results are 
reported in and will be described more fully in |14j . 

Appendix A: other rotation centres 

In the main text, we have analyzed the standard situation of coincidence rotations 
around lattice (or module) points. Here, we will briefly comment on rotations around 
other centres in the lattice case and on situations with more than one translation 
class of points. 

n=4: the square lattice revisited 

Another obvious rotation problem is that around the centre of a Delaunay cell of Z^, 
(|, |) say. This point represents the only class of deep holes of Z^, cf. [Ej, and has 
the entire point group of Tl? as site symmetry. It is obvious that the coincidence 
problem is equivalent to that of the point set F defined by 

r = {a + i6|a, 6gZ, a + 6 odd} , (Al) 

which is obtained from Z^ — (|, |) via rotation through 7r/4 and dilation by 
Observe that ()A1|) can be rewritten as 

r = {« G Z[i] I « ^ (l+i)} (A2) 

which solves the problem: as was shown in Section 2, the coincidence rotations of 
T? can be factorized, the generators being e^"^ = i (rotation through vr/2) or of the 
form e*''' = ujp/uJp with N{ujp) = p = 1 (4), hence Up ^ (1 + i). The former still 
is a symmetry of F (index 1), and we also get the latter because both numerator 
and denominator are in F. Also, the reflection in the x-axis remains a coincidence 
operation of index 1. Summarizing: 

OC(r) = OCiT?) , (A3) 

and the coincidence indices are unchanged. 

n=3: the hexagonal packing 

Consider the Voronoi complex of the triangular lattice — it is a packing made from 
regular hexagons — and let H be its vertex set. Let us consider rotations around 
the centre of a hexagon which is a point of maximal site symmetry Dg. If we rotate 
the complex through 7r/6 and dilate by -\/3, then R can be characterized as 

E = {a(i %\q\ I a ^ (1 + ^)} (A4) 

where ^ = i+iv^. Since A^(l + f)) = 3 and 3 is not a complex splitting prime in Q(^), 
we find again all rotations and reflections which we had already for the triangular 
lattice: 

OC{R) = 0C{A2) , (A5) 
and also the indices remain unchanged. 
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n=3: coincidence definition revisited 

Slightly different is the situation if we keep the entire set of lattice points, but rotate 
around the centre of a Delaunay cell: the latter is a triangle and its centre has only 
site symmetry. We rotate again through 7r/6 and dilate by which gives the 
point set 

G = {ae %\q\ \ (x = \ (1 + ^)} . (A6) 

Here, a rotation through 7r/3 would change the congruence class of G from 1 to — 1, 
so it is no longer a coincidence rotation. This reduces the torsion part of OG from 
Cg to C3 in agreement with the reduced site symmetry while all other generators 
remain unchanged. In particular, the reflection in the x-axis leaves G invariant and 
the index formula applies for all remaining elements. 

One might also consider possible variants of the coincidence concept here: a 
rotation through 7r/3 alone does not produce a coincidence for the set G, while the 
same rotation followed by a suitable translation can give a coincidence of index 1. 
The latter might be more important when the connection to grain boundary growth 
is considered. Indeed, especially in view of applications to nonperiodic discrete point 
sets, one might define (with obvious meaning) 

inf[P : Pn(i?P + t)] (A7) 

to be the coincidence index of an isometry R acting on a point set P. This gets rid 
of the dependence of the index on the rotation centre and comes closer to the idea 
of optimal fitting of grain fragments. 



n=5: the rhombic Penrose tiling 

A complication here is that the vertex sites of the rhombic Penrose tiling T fall into 
4 different translation classes with respect to the uniquely defined limit translation 
module A1(T), compare (23 1201 ■ identify A1(T) with the projection of the 4D 
root lattice into tiling space for definiteness. Then each point class has its own 
window of pentagonal shape. The windows come in two different sizes (related by a 
factor of r = (1 + -\/5)/2) and in pairs related by rotation through vr, compare |29] . 
The vertices of the rhombi are not points of the module 7V1(T) (which also means 
that none of them is a "standard" rotation centre). 

Let us now consider the coincidence problem of the set of vertex sites with all 
translation classes identified. To be explicit, we take the rhombic version of the 
cartwheel pattern where the rotation centre is not a rhombus vertex but coincides 
with the centre of a regular decagon filled with rhombi. This point is a representative 
of the fifth translation class, so far absent. The cartwheel tiling has Dio symmetry 
in the sense that any Dio-operation either maps the tiling upon itself (thus, in 
particular, the set of vertex sites) or produces at most a mismatch of density zero 
(along the well-known worms). All these operations thus have coincidence ratio 1. 
The corresponding rotation in window space maps windows to windows, because 
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they appear in Dio-orbits around the origin. More than this, it maps translation 
classes of windows to translation classes of windows. 

For other coincidence isometrics, we first observe that the the integral span of all 
vertex points is again a planar module of rank 4, in our explicit case the projection 
of the weight lattice Al, the dual of A^, into tiling space. This module is equivalent 
to M.{T) and possesses therefore the same coincidence isometrics, namely those 
described in Section 3. Consequently, we find all these also as coincidence isometrics 
of the rhombic cartwheel tiling. The coincidence ratio must now be corrected in a 
similar way to that of the Ammann-Beenker tiling in Section 4, but the window 
system requires a slightly more complicated calculation which we will not present 
here. 

Even more complicated would be the coincidence analysis for rotations around 
vertex points, in particular with various point classes distinguished. The methods 
needed are in principle those described forn = 3 above, but details will not be given 
here. 

n=12: a square-triangle tiling 

Quasiperiodic square-triangle tilings are attractive for a number of reasons. We 
mention them because they can have 12-fold symmetry in the sense of mismatches 
of at most density zero under Di2-operations, see for an example. There, all 
vertex points are in one translation class, so no problem occurs and we find all 
coincidence isometrics of Section 3. But for the correction factor due to window 
overlaps one encounters a new type of complication: the window is fractally shaped 
and consequently we see no way of calculating this factor. It is left as an exercise 
for fractal readers. 

Appendix B: proofs 

Here we give the promised references and proofs of Facts 1-3 in Section |21 

Fact 1 This is proved in (2^ (Lemma 4) or j2HI (Theorem 2.13) for example, 

but we sketch the proof here as it leads on naturally to the proof of Fact 2, which 

is less commonly found in the literature. 

Let P be a prime factor of p in K. Since p\n the nth roots of 1 in are distinct 

mod P. (The most straightforward way to see this is from the identity 

n— 1 

n = X{{l-e), (Bl) 

k=l 

got by putting a; = 1 in (x" — l)/(a; — 1), and noting that every difference of roots of 
unity is an associate of 1 — for some k.) The residue class field Fp = Z[^]/P is a 
finite field generated over by the residue class ^* of ^, and since distinct roots of 
unity are distinct mod P, the order of ^* in Fp is n. Every finite extension of F^ is 
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normal with cyclic Galois group generated by the Frobenius automorphism x t-^ x^, 
whose order is the degree d of the extension. Consequently the degree [¥p : F^] 
is the smallest d with = that is, the smallest d with n \ {p'^ — 1). This 
establishes Fact 1, since \¥ p : F^] is the degree of the minimal polynomial satisfied 
by ^ mod p and hence is deg^(p). 

Fact 2 Analogously to the above proof, degj^ip) is the degree d' of the residue 
class field extension F^/F^, where p is the prime of L divisible by P and F^ is the 
residue class field + ^^P- Clearly [Fp : Fp] <[K : L\ = 2, so the d of Fact 1 
is either d' or 2d' . If d is odd then d' = d, the order of p mod n, and no power of p 
is congruent to —1 mod n. 

To treat the case of even d we first note that if C,*''^ + = ^*'=2 + ^*-fe2 jp^ 
(where < ki, k2 < n) then either ^*^^ = ^*''^ or ^*^^ = ^*-^^. This is because ^*''^ , 
Qj-e the two roots in x of 

^2 _ ^ ^ ^ ^ ^ 2) (B2) 

and when the equations are the same the roots must match in some order. Now 
d' = d/2 if and only if ¥p is the unique sub field of index 2 in Fp, this being the 
fixed field of the element x i— > x^''^^ of order 2 in the Galois group of Fp/F^. So 
d' = d/2 if and only if 

^* + f -1 = + = ^*p"/' + C~^'^\ (B3) 

which requires ^*~^ = ^*p'^^^ (equivalent to n \ p'^^'^ + 1), since ^* 7^ ^*p''/^^ 'jj^g 
exponent (i/2 here is plainly minimal, since n\p°' + l^n \ p^"" — 1. 

Fact 3 Part (a) is a result of the fact that when n = p^ then p is totally ramified 
in K (that is, p is the 0(p'')th power of a degree 1 prime of K), see for example |23], 
Lemma 3. As a consequence deg^(p) = 1 and p is not a complex splitting prime 
because it has only one prime factor in K. 

For part (b) we refer to the Hasse diagram of field inclusions in Figure 4. Here 
Ki and K2 are the cyclotomic fields of riith and p^th roots of unity and Li and L2 
their maximal real subfields. Then K = K1K2, the compositum of Ki and K2, and, 
since p /ni, fsTi n i^a = Q (see Thm. 9.52 or [28j Prop. 2.4). Let p = be a 
prime of K dividing p. For an arbitrary subfield F of K we denote by pp the prime 
ideal of F that is divisible by p. 

Because p is unramified in Li and Ki but totally ramified in K2 it follows that 
and are totally ramified in K2L1 and K and, in particular, 

degi^,i,/i,(PK2Lj = deg^/^^(p^) = 1. (B4) 

Consequently 

degK/K^L^ (P) = deg^^/i^ {Pk^ ) . (B5) 

Now look at the fields L1L2, K2L1, L and K. Since p^^ ramifies in K2 but p is 
unramified in Li, p^^i^ ramifies in K2L1 and hence in K. By Prop. 2.15(b) of ^H] . 
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K 




Pi is unramified in K, and hence Pl^l^ ramifies in L. We now have 

degL/L^L-APL) = degK2Li/LiL2(Pi^2Li) = 1 (B6) 

whence 

degi^/L (p) = degK/K^L, (p) • (B7) 

Since Pj^ is unramified in K and Pj^^ is unramified in i^'i, Eqs. ()B5|) and ()B7|) imply 
that Pj^ factors into two primes of K if and only if p^^ factors into two primes of 
Ki. Finally, deg^(p) = deg;^^(j9) is an immediate consequence of the fact that pj^_^ 
is totally ramified in K. 
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